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SUMMARY 


The flow about a body of revolution at hi^ supersonic airspeeds is 
investigated analytically with the aid of the generalized shock-expansion 
method. With the assumption that flow at the vertex is conical, approxi- 
mate solutions for the flow field are obtained for values of the hyper- 
sonic similarity parameter (i.e. , the ratio of the free-stream Mach number 
to the fineness ratio of the body) greater than about 1 and for angles of 
attack less than the semivertex angle of the body. Surface streamlines 
are approximated by meridian lines and the flow field is calculated 
in meridian planes. Simple explicit expressions are obtained for the 
surface Mach numbers and pressures in the special case of slender bodies. 

In the case of lifting cones , algebraic solutions defining the entire 
flow field are obtained when the hypersonic similarity parameter has a 
value of about 1.4 or greater. 

Surface pressures and shock-wave shapes were obtained experimentally 
at Mach numbers from 3*00 to 5*05 and angles of attack up to 15^ for two 
ogives having fineness ratios of 3 and 5 and for two cones having the 
same vertex angles as the ogives. The predictions of the methods of this 
paper are found to be in good agreement with experiment at values of the 
hypersonic similarity parameter in the nei^borhood of 1 and greater, 
when the ratio of angle of attack to semivertex angle is about one-ha] f, 
or less. For increasing values of this ratio, agreement deteriorates but 
may still be considered fair for values sli^tly less than 1. 


INTRODUCTION 


It was suggested in reference 1 that flow over the surface of a non- 
lifting body of revolution could be treated as two-dimensional in type 
downstream of the vertex when the hypersonic similarity parameter (i.e., 
the ratio of the free-stream Mach number to the fineness ratio of the body) 
was greater than about 1. This point was substantiated by comparing predic- 
tions of two-dimensional (Prandtl -Meyer) expansion theory with those of 
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characteristics theory for the Mach numbers and pressures on the surfaces 
of ogives. The tvo -dimensional theory has the advantage^ of course^ of 
being relatively simple by comparison to characteristics theory and is 
about as simple as the recently proposed hypersonic small disturbance 
theory of Van Dyke (ref. 2 ) . 

It vas also suggested in reference 1 that the tvo -dimensional approach 
might be extended to the calculation of flow at the surface of sli^tly 
inclined bodies of revolution. This thought led to a study (ref. 3) of 
three-dimensional hypersonic flows which revealed that such flows may 
often appear locally two-dimensional. It was concluded that at hypersonic 
speeds the entire flow field about a three-dimensional body may, under 
certain conditions, be calculated with a shock-expansion method similar 
to that employed for calculating two-dimensional flow about airfoils 
(ref. h) , The conditions of when and how this generalized shock-expansion 
method can be applied to calculate three-dimensional flows were determined 
in reference 3 * 

The principal objective of the present paper is to apply the general- 
ized shock-expansion method to obtain expressions yielding the Mach nximber 
and pressure distributions throu^out the entire flow field about an 
inclined body of revolution. In order to apply the shock-expansion method, 
it is necessary to know initial conditions at the vertex of a lifting 
body. These conditions can be taken to be the same as those about a cone 
tangent to the body at the vertex. One objective of this paper, then, is 
to develop a conical flow theory for lifting cones over the range of free- 
stream Mach numbers and apex angles not treated in the M.I.T. tables 
(ref. 5 ). 


NOTATION 

a local speed of sound, ft/sec 

^ . normal force 

Cv normal-force coefficient, 

Cp specific heat at constant pressure, ft-lb/slug ^ 

Cy specific heat at constant volume, ft-lb/slug 

d maximum diameter of body of revolution, in. 

E entropy, ft-lb/slug ^ 

H total pressure, Ib/sq in. 

K hypersonic similarity parameter, Mq y 

I characteristic body length (measured from vertex to most for- 

ward point of maximum diameter) , in. 
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Mach number (ratio of local velocity to local speed of sound) 
static pressure^ Ib/sq in. 

dynamic pressure^ Ib/sq in. 

P - Po 

pressure coefficient, — — — 
gas constant, ft-lb/slug ^ 

velocity component parallel to ray passing throu^ vertex of 
cone, ft/sec 

velocity component normal to u in a meridian plane, ft/sec 
velocity component normal to a meridian plane, ft/sec 
resultant velocity, >/u^ + v^ 4- , ft/sec 

maximum velocity obtainable by expanding to zero temperature, 
ft/ sec 

distance along axis of body measured from vertex, in. 

distance normal to axis of body, in. 

angle of attack, radians unless otherwise specified 

Mach angle, arc sin ^ , radians unless otherwise specified 

C-p 

ratio of specific heats, 

Cv 

angle of flow inclination in meridian plane measured with respect 
to body axis, radians unless otherwise specified 

angle of inclination of axis of conical shock with respect to 
free-stream direction, radians unless otherwise specified 

angle of inclination of axis of conical shock with respect to 
axis of body, radians unless otherwise specified 

mass density, slugs/cu ft 

angle of merdian plane with respect to plane of symmetry, 
radians unless otherwise specified (see fig. l) 

angle between axis of cqne and ray passing throu^ vertex of cone, 
radians unless otherwise specified 
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Subscripts 


o 


free-stream conditions 



conditions at different points in the flow field 


c 


conditions on the surface of a cone 


e 


conditions at the external surface of the vortical layer 


conditions on the surface at the vertex of a body 


s 


conditions inmediately behind the shock wave at the vertex of a 
body 


THEORY 


This investigation is concerned with the theoretical and experimental 
characteristics of the flow about bodies of revolution traveling at high 
supersonic airspeeds and at small angles of attack, it is assumed through- 
out the analysis that the disturbed flow is everywhere supersonic and, 
hence, the body has a pointed nose or vertex. With these restrictions 
on the free-stream Mach number, angle of attack, and body shape, the bow 
shock wave will lie close to the surface of the body. The procedure for 
determining flow conditions in such flow fields is analogous to that 
employed in reference 1 for the case of axially symmetric flow fields; 
namely, the flow field is studied in two parts - flow at the vertex and 
flow downstream of the vertex. The combined results of these two phases 
of the investigation will then be applied to the determination of the 
whole flow field and, in particular, to the determination of flow prop- 
erties on the body surface and the resulting shock-wave shape. 


It follows from the assumptions basic to this analysis that the flow 
at the vertex will be the same as for a cone tangent to the body at the 
vertex and immersed in the same free stream. All derivatives with respect 
to radial distance vanish for these conditions, and the equations of motion 
and continuity in spherical polar coordinates become (a schematic diagram 
of the polar coordinate system is shown in fig. l) 


Flow at the Vertex of a Lifting Body 
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bw ^ V dw 
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1 
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(la) 


(lb) 


(Ic) 


and 


2pu sin (0 + V 


dp . dv dp 

sin to + p sin (0 — + vp cos to + w — + 

dto dto dtp 


p|s = 0 

dep 


( 2 ) 


respectively. Since the total energy in the flow is constant, the fol- 
lowing relations must be satisfied: 


r A ^ 

7 - 1 \P dtp 


depy 





•y - 1 \P dto p^ dto/ \ dto dto dto / 


(3) 


The entropy at any point in the flow may be expressed as 


^ In 

fp 

r - 1 

1 

An 

a 

1 


Equations (l), (2), and ( 3 ) together with the relation 


(i^) 




V^) 


may be combined (by eliminating the pressure and density terms) to obtain 
the general equation of motion 
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uv2 - uw2 


by __ by _ w 

Su) sin CO Sep sin w Sep 



which, when combined with the appropriate shock-wave relations, defines 
the flow about a circular cone immersed at an angle of attack in a super- 
sonic stream. 

Conditions at the shock wave .- In order to obtain algebraic solutions 
for flow at the vertex, it is necessary to make some simplifying assump- 
tions regarding the flow field. To this end, the conical shock^is assumed 
to be circular but inclined at an angle e to the free stream. ^ Then, 
the angular difference between the cone axis and the shock axis is 


Now, the shock-wave angle measured from the cone axis, cos, is referred 
to the free-stream direction by cos + ex cos cp. The shock-wave angle 
measured from the shock axis, (ws)cp=o referred to the free-stream 

direction by (u)s)cp=o + “H + ^ ^os cp. Hence, 


and by virtue of expression (6) the resulting equation for the shock 
angle may be written^ 


where (us)cp=o is the shock angle in the plane of symmetry on the windward 

^Experimental results indicate (as will be shown later) that for small 
angles of attack the conical shock does, in fact, remain nearly circular. 
Other investigators (notably Stone, see, e.g., ref. 6) have made a similar 
assumption. It should be noted, however, that the additional assumption 
commonly employed, namely, that the conical-shock apex angle is the same 
as in the nonyaw case, is not made here. 

^It should he noted that ail angles are measured with respect to the 
body axis. The procedure of developing all pertinent expressions in the 
body coordinate system will be employed throu^out the analysis. 


q = a - € 


( 6 ) 


u)s + a cos cp = (ws)cp=o + T + ^ cos cp 


= (<Js)cp=o + - cos cp) 


( 7 ) 
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side of the cone (see fig. l) . Similarly, the meridian angle cpg, measured 
with respect to the shock axis, may be related to the meridian angle cp, 
measured with respect to the cone axis, by 


sin cpg = sin cp(l - t] cot o) cos cp) 


and 


cos cpg = cos cp + T) cot 0) sin^ 

The shape of the conical shock having been specified, the velocity com- 
ponents at the shock may be obtained from the oblique-shock-wave relations 


Ug = Vocos[t0g + a(cos cp + T] cot cogsin^) ] (8) 

= 7 q€ sin cp(l - T cot cjgcos cp) (9) 



/N 

V 


7-1 

7 + 1 ^ 

V 



sin[tjg + a( cos cp + t) cot cjgsin^) ] 


( 10 ) 


If the shock angles in the plane of symmetry were known, flow conditions 
around the entire shock could be determined since the shock-wave angle 
in any plane could be determined from equation (7), In order to determine 
these shock-wave angles , it is necessary to determine the crossflow 
component of velocity, w, throughout the flow field. Attention is there- 
fore turned in this direction. 

Determination of crossflow component of velocity. - Recalling the 
basic assumption to this analysis, namely, that flow fields of the type 
under consideration are characterized by the bow shock lying close to 
the surface of the body (i.e., oj -6 <<1), it is reasoned that the 
variation of w with u should be small and it is assumed that 


V = wi(cu)w2(cp) (11) 

Now ws(cp) is given by equation (9). There remains, then, the determina- 
tion of wi(w) in any meridian plane. To this end, consider equation (la). 
Differentiating this equation twice with respect to cp and once with 
respect to u yields, in the plane of symmetry. 
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Sf') . sin u ('Ai + ^ cot u> 
/cp=o \ScpSw 


( 12 ) 


Now, if we let 


^ '' = F(») 




(13) 


equation (12) becomes, upon integration with respect to w. 






— sin CO = ( ) sin cOg - / F((o)dco 




o 


0) 


(14) 


Consider now the integral term in the above equation. At the surface of 
the body (v = O) 


F(co) =0 

Since it seems reasonable that F(io) will be a monotonic function between 
the surface and the shock wave, it is assumed that this fraction will 
attain its maximum value at the shock and may be written (from eqs. (o) , 
( 9 ) , ( 10 ) , and ( 13 )) 


F(u))s = 


2(y - 1 ) 

7 + 1 


Vocos(cos +a-)(e - 2ox\ cot (Og) + 


(Vo6)Mo€ 
Mosin(c0g + o,)_ 


7-1 

7+1 


1 + 


7-1 ,, 2 „. 2 


Mo sin^(cos + 


(e - 2aq cot cOg)VoCOs(tOg + ct) 


Now, according to the above expression, F(co) is a maximum when 
M sin((0g + oc) is a minimum. Since Mosin(tos + a) > 1 for attached shock 
waves, the maximum value of F(co) can be determined by setting 
Mosin(u)s + cc) equal to 1 in the above expression. Hence, since Mo€ ~ 1 
there results 


F(oj) 


2(7 - 1 ) 


Vo€ 


3 ~ . 1 Vocos(o)s + a)(€ - 2ot,Ti cot Ug) 
7+1 


7+1 
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Hence, 

|f(w) I < e 

and it follows that 


-^U)c 


F(u))da) < 


.Cl)s 


e dco = €(cjs " Wc) 


Cl) 


CJr 


Therefore the integral term in equation (lU) is, at most, second order 
and can he neglected in this analysis. Equation (l 4 ) now can he written 


^w _ ( ^ 

Sep \Scp/g sin Cl) 


Hence ^ 


Wi(w) 


sin u)s 
sin w 


and it appears that expression (ll) is a logical assumption since the 
variation of w with w is, in fact, small. Combination of the above 
expression with equations (9) and (ll) yields, then. 


sin (jjg , , 

w = Yq£ ■ sin cp(l - T] cot Wgcos cp) 


(15) 


defining the crossflow component of the velocity anywhere in the flow 
field relative to the cone axis. 

Having determined an expression for w throu^out the flow field, 
one can now obtain a solution to the flow in the plane of symmetry in a 
manner analogous to that presented in reference 1 for the case of axially 
symmetric conical flows. Since the calculation of the flow in this plane 
requires simultaneous solutions of the conical flow equations and the 
oblique shock-wave relations, the procedure is somewhat involved and, 
hence, is given in Appendix A. After (u)g)q)_Q and € have been determined 
from Appendix A, conditions around the entire shock front can be deter- 
mined from equations ( 6 ) throu^ (lO). Determination of flow conditions 
around the cone surface will now be considered. 

Flow conditions on the surface. - It has been shown by Ferri in refer- 
ence 7 that to the first order in a the entropy remains constant in a 
meridian plane (having the value that exists at the shock in that plane) 
until a vortical layer is reached at the surface of the body across which 
a variation of entropy occurs. Since the entropy on the surface is 


10 


MCA TN 33U9 


constant, it must have the same value that exists in the plane cp = 0. 
Now the thickness of the vortical layer being of the order of a^, the 
pressure remains constant through the layer to the first order of a 
and the component of velocity normal to the surface can he considered 
zero on the external side of the layer. Hence, across this vortical 
layer 


Pc = Pe 


V 


c 


= 0 


and 


Ec - (E)(p=o 


where the subscripts c and e refer to quantities inside and outside the 
layer, respectively. Consider now an expression relating u and w on 
the surface of the cone which may be obtained from equation (la); namely. 


— = w sin 5. 

acp 


or 


ncp 

^ = (Vc)cp=o + / ’W ‘ 3 .CP (16) 

o 

Since the thickness of the vortical layer is proportional to a^, the 
normal component of the velocity is zero through the layer and the above 
expression holds on either side of the vortical layer. An expression for 
the velocity and, hence, the Mach number externally adjacent to the vorti- 
cal layer may then be easily determined as follows. The expression for 
the crossflow at the vortical layer may be written in the form (from 

eq. (15)) 




sin u)g 

sin cp(l - T] cot cjgcos cp) 


(Vc)cp=o Sc 

Substituting this expression in equation (l6) yields^ upon integration^ 


(IT) 


u, 


(Vc)cp=o 


- =1. 


(Vc)cp=o 


Sin o)c 


cos cp - cot cjgsin^cp ) (l 8 ) 


e 
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The Mach number at the external side of the vortical layer. Me, may be 
determined from the relation 


Me^ = 


(Mc)J=o 


L(Vc)(p=oJ 




f 

1 

(L) 

> 

1 

/ 

r 

1 

CV) 

1 

o 

II 

o 



(19) 


where, of course. 



2 2 
ue + Ve 

(Vc)^=o 


The velocity and the Mach number in the cp = 0 plane are given by ( see 
Appendix A) 



^s“^s 


cos(cl)s - Sg) 

(be - b) (o)s - bg) 


cos 


and 


(M) 


2 

cp=o 



respectively, where 


(7 + l)^Mo‘^sin^ 
(^s)ip=o ~ r 

4 2/Mo' 


If one employs the condition that the pressure is constant across the vortical layer, an 
expression defining the Mach number directly on the surface of the body. Mg, in terms of the Mach 
number at the vortical layer, may be obtained; namely. 


(‘^s)cp=o 


- 4 -j^MQ^sin^ 


(‘*^s)cp=o ^ 


l^ ^/Mo^sin^ 


(‘^s)cp=o + ^ 


^ V 


ro 


■sin^ 


(wc,)^ „ + a 
s Cp=o 




(r - 1) H " l)Mn^sin^ 


(wc.) + a 

^ ® cp=o 


.2} 


yir 


7 - 1 ^ 


^e”^c 
+ Mo^le 4'Cv 


7-1 


(20) 


where (since Ec = (E) 0 =q) 

Eg-Ec 

7Cy 


4 o^sin^ 

+ a 

1 

1^(7 - l)Mo^sin^ 

o)s 

-1- a(cos cp -H T] cot cogSin^cp) 

-I 

\^^sin^ 

u)g -t- a(cos cp - 

_ 

f T] cot (A)gSin^) 


j ^(7 - DMo^sin^ 

(tOo) + a 

s cp=o 

^ 4 , 


'7 


27Mo^sin^[o)s + a(cos cp + r) cot cogsin^cp) ] - (7 - l) 


27Mo^sin^ 


( U)q ) OL 

^ ^^cp=o 


(7 - 1 ) 


(21) 


Since the flow is isentropic in the cp = 0 plane (downstream of the shock wave) and around the 
surface of the body inside the vortical layer, the pressure coefficient anywhere on the surface 
may, of course, be obtained by the expression 
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op=o 


1 


( 22 ) 


where 



(p=o 


7 + 1 


(7 - 1) 






and is given by equation ( 20 ) 


Floy conditions off the surface . - Flow conditions in the plane of 
symmetry, on the surface of the cone, and at the conical shock having 
been determined, the flow throughout the remainder of the flow field may 
now be calculated. Since only hi^ Mach number flows are considered in 
this analysis, the variation of the magnitude of the resultant velocity 
in a meridian plane will be small. Hence, the variation of u and v will 
be small and may be represented by a power series in (u - &c) where the 
coefficients are determined by the requirements 


u = Ug at CJ = 5 (, 


u = Ug at u = lOg 




and 


V = 0 at 0 ) = 


at w = tUg 




Ik 
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Thus , 



where Wg, us, vs, and ue are given by equations (?), (8), 

(18), respectively. There remains now the determination of (du/d 
(8v/8w) g. To this end, consider equation (lb) . Just outside the 
layer this equation reduces to 


= we^cot 6g 

It will be recalled that the entropy was assumed constant between 
shock and the vortical layer in each meridian plane. Now Euler’s 
tion for compressible flow along a stream tube may be written 
(cp = constant) 



1 ^ = -V ^ 
P 


Combination of this expression with equation (25) yields 


(v 1^^ = - Wg^cot 5g 


from which may be obtained (noting that at the surface 


= u^ 



where Wg is given by equation (l?) and 


8wN _ _ V e cot 6gsin cp(l - t| cot Wgcos cp) 

.8o)/g ° sin 6c 


(23) 


(24) 

and 

a) g and 
vortical 


(25) 

the 

equa- 

w^) 

(26) 

(27) 
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Wow from equation (5) there results (setting v = O) 




Ue + CSC m 


Me' 


2/We 


- 1 - u. 


(28) 


where Me is known from equation (19) and from equation (I5) 


= Voe (cos cp - q cot cjscos 2 cp) 


sin u) 


(29) 


The components of the local velocity anywhere in the flow field external 
to the vortical layer are now known from equations (15) j (23), and ( 24 ). 
Hence, the magnitude and direction of the resultant velocity and, con- 
sequently, the Mach number may easily be determined. If the Mach number 
is known, the local pressure (in coefficient form) may be obtained any- 
where in the flow field with the aid of the expression 


P 


2 _ /£s _p_ 

7Mo^ V^° Ps 



(30) 


where 


h. 

Po 


27Mo^sin^[a)g + a(cos cp + q cot Ugsin^cp) ] - (7 - l) 

7+1 


and, of course. 



(31) 


(32) 


The Mach number and pressure distribution (as well as the orientation of 
the conical shock) are now known throughout the flow field around a lift- 
ing circular cone. 

The range of applicability of the results of this analysis is con- 
sidered to be the same as that of the nonlifting cone solutions presented 
in reference 1 . This results from the fact that when a = 0 , equation (AT) 
in Appendix A reduces identically to the equation defining the deflection 
angle in reference 1 , As was pointed out in this reference, when 
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Mg(us - 5c) > 1/2 an Imaginary value of 5s is obtained. If a > 0, 
equation (A 7 ) yields a real value of 6 g, (for Ms(wg - 5c) > I/ 2 ) , how- 
ever, it would not be expected that the range of applicability of this 
equation (in terms of Mq and 5c) would be increased for finite a. 

Figiire 2 shows the boundary line (given by Ms(Wg - 5c) = l/2 for a = 0 
above which the present conical flow solutions are applicable. The dashed 
line represents the boundary below which the results of Stone's second- 
order solution (ref. 5 ) are available. 

The flow around circular cones traveling at small angles of attack 
and at high supersonic airspeeds can be calculated by means of the fore- 
going algebraic expressions. As was pointed out previously, these expres- 
sions can be employed to determine fluid properties at the vertices of 
pointed bodies of revolution other than cones. Investigation of flow 
downstream of the vertices of such bodies is now undertaken. 


Flow Downstream of the Vertex 


In this study we exploit the finding of reference 3 that many three- 
dimensional hypersonic flows may be treated by a generalized shock- 
expansion method which is analogous to that employed in reference 4 for 
two-dimensional flows. Specifically, this treatment is permissible when 
disturbances associated with the divergence of streamlines in planes 
tangent to a surface can be considered negligible compared to those 
associated with the curvature of streamlines in planes normal to the sur- 
face. For the case of noninclined bodies of revolution which are curved 
in the stream direction, this requirement is satisfied when the hyper- 
sonic similarity parameter is greater than about 1 (see ref. l) . For 
inclined bodies, an additional restriction is imposed. This point is per- 
haps best clarified by considering the problem of calculating flow at the 
surface. 

It follows from reference 3 that when the generalized shock-expansion 
method applies in the region downstream of the vertex, surface streamlines 
can be approximated by geodesic lines. The only geodesics on the sxirface 
of a body of revolution which, like streamlines, do not intersect each 
other are the meridian lines. In addition, the meridian lines are the 
only geodesics which, like the streamlines, pass through the vertex. When 
the shock-expansion method is applied, then, surface streamlines are 
approximated by meridian lines. Strictly speaking, however, this can be 
the case only when a < < 1. (It should be noted that this is always 
true, independent of <x, for the extreme windward and leeward streamlines.) 
Evidently, then, the generalized shock-expansion method should be appli- 
cable to curved bodies of revolution only at small angles of attack in 
flows characterized by a value of the hypersonic similarity parameter 
greater than about 1. 
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The procedure for determining flow conditions at the surface of a 
lifting body is entirely analogous to that employed in the application 
of the shock-expansion method to the nonlifting body (ref. l) . For 
example^ the Mach number on the surface at the vertex is obtained with 
the aid of equation (20) . The variation of Mach number downstream of 
the vertex is then obtained by means of the isentropic expansion relation 
(see ref. l) 


8b - Sa = sin"^ i 

Ma 


Sin 


-1 


M- 


■B 




J 


7 + 1 


tan 


-1 


7+1 < 

(7 - 1)(Mb2 - 1) 



7 + 1 

- 1)(Ma^ - 1) 


( 33 ) 


where A and B are different points on the same meridian line (or stream- 
line) . If the Mach number distribution is known^ the pressure distribu- 
tion (in coefficient form) on the surface is readily obtained with the 
aid of equation (22) . It should be noticed that when Me is employed^ 
equation ( 33 ) yields the Mach number distribution on the body under the 
vortical layer This result materially reduces the net labor associated 
with carrying out the calculations to determine the pressure distributions 
around the body downstream of the vertex since the pressure rise across 

the shock need be considered only in the plane cp = 0 (see eq. (22)). It 

dCj^ 

results, too, that a relatively simple expression for - 5 — , the initial 

a=o 

slope of the normal-force-coefficient curve, can be obtained. The devel- 
opment of this expression follows. 


The expression for the normal-force coefficient on a body of revolu- 
tion may be written 


Cn 


16 



Po 


r cos cp dep dx 


(3h) 


where d is the diameter of the base and r is the local radius of the 
body. Differentiating expression (3^) with respect to a and employing 
the condition of constant entropy on the surface results in 


^It is reasoned that since a vortical layer exists around the body 
surface at the vertex^ then a vortical layer must exist downstream of 
the vertex as well. 
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dCN 


da 


a=o 


16 

rMo^ifd^ 


J J W/ V o^ 

o o 


r cos cp dp dx 


(35) 


This expression may be rewritten in terms of the initial normal-force 
curve slope for a cone tangent to the body at the vertex; namely, 


dCN 


da 


a=o 


= 8 


dj da 


tan 6m 

TCN 



/sin 2Pu\ 

\sin 2P J 



(36) 


where the subscript TCW refers to a cone tangent at the vertex. The 
calculations necessary to determine the initial normal-force-curve slope 

for a body of revolution are relatively simple, since may easily 

be obtained from Stone’s first-order theory (ref. 6) or from chart 8 in 
reference 8. The Mach number and pressure distribution along the body 
may be obtained by the conical shock-expansion method presented in refer- 
ence l.'^ Having determined these distributions, one may easily evaluate 
the integral term in equation (36) by numerical integration or by graphical 
methods . 


In order to determine fluid properties in the flow field other than 
on the surface, it is necessary to know flow conditions just outside the 
vortical layer downstream of the vertex. These conditions may be deter- 
mined in the same manner as before except that now initial flow condi- 
tions externally adjacent to the vortical layer at the vertex are 
employed. For example, the Mach number at the vertex is determined by 
means of equation (I9) . Equation (33) is then employed, as before, to 
obtain the Mach number distribution downstream of the vertex. When flow 
conditions along this layer in a meridian plane have been established, 
fluid ‘properties throughout the plane may be calculated after the manner 
described in reference 4. Application of the generalized-shock-expansion 
method for determing the flow field in any meridian plane is discussed 
in Appendix B. 


Simplified Expressions for Slender Lifting Bodies 


In the case of slender bodies traveling at very hi^ Mach numbers 
(again a < < l) the calculations of fluid properties at the surface 
become relatively simple and, hence, merit special attention. As in the 
case of the nonlifting body, a hypersonic slender-body theory yielding 

^It is clear that the shock-expansion method discussed previously in 
the present paper may also be used since the expressions developed herein 
reduce identically to those of reference 1 when a = 0. 
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explicit solutions for the Mach number and pressure at any point on the 
surface of the body can be obtained. These solutions may be sxmimarized 
(from Appendix C) as follows. The local Mach number at any point is given 
by the relation 


M = 


% 

1 - 2-^ (MnSb)(^i - 


(37) 


where 






(Ms5]\j) 


1 . 


Wc 


^“2^ (%e)p=c 


2 (jJc 



and 



n 2 


^ cp=o 


(7 + 1)% (<^s)(p=o + 


27 M 0 


n 2 


(“s)cp=o + "" 


- (r - 1)}|(7 - l)Mo^ (ws)cp=o + 2 } 


Ec ■ Ee 
7Cv 


-,2 


2rMo (ws)(p=o 


- (7 - 1 ) 


1 

7 


27 Mq^(cjs + a cos cp)^ - (7 - l) 


Mo^(o)s + O' cos cp) ^(7 - l)Mo 



-i 2 


(u)o) + a 


+ 2 


n 2 r 


(“s)cp=o + 


2 2 
(7 - l)Mo (cOs + a cos cp) +2 


Unless otherwise designated^ 0)3 in the above expressions is given by 


<^s = (“s)q,=o 


+ T](l - COS Cp) 
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and (ws)q, 3 ,o determined from equation (CT) in Appendix C. Now 


7 + 1 ^ 
4 &N 



D 


a 


0 


where 



and, of course 


€ = a - T| 


The pressure coefficient at any point on the surface may he obtained from 
the expression 



(^s^cp 0 _ 7 - 1 5 ) ^ 

Mu 2 ® ^ T-oV 

27 

- 1 

^ ^cp=o 




It is interesting to note that equations (3?) and ( 38 ) predict the ratios 
of local to free-stream Mach numbers, and local to free-stream static 
pressures to be the same at corresponding points on related bodies, pro- 
vided that the flow fields about these bodies are related by the same 
respective values of the hypersonic similarity parameters Mo5jf and Moa. 
These predictions® are identical to those of reference 9 for inviscid 
flow about slender three-dimensional shapes. Hence, these expressions 
readily lend themselves to solutions in terms of Mq5n 9^*1 tab- 

ular form. Calculations over a range of Mq 6 n from 0.6 to 3 .O and a/Sm 
from 0 to 1 were carried out for flow at the vertex of a body of revolu- 
tion and the results of these calculations for the flow parameters 


increments of cp from 0° to l80°. Thus, for a given Mo 6 u and Moa, the 
Mach number on the surface of a body downstream of the vertex is readily 
obtained with the aid of these tabulated parameters when used in conjunc- 
tion with equation ( 37 ). The pressure coefficient anywhere on the surface 

of the body is e asily calculated by means of equation ( 38 ). 

®In the case where a = 0 , the expressions developed in the present 
paper reduce identically to those presented in reference 1 . 


(Ms)cp=o (Ps/Po) 
% ^ (Mo5w 


, and (MsSn)|^_q are tabulated in table I for 30° 


MCA TN 33h9 


21 


The results from table I may also be used to good advantage in deter- 
mining the initial normal-force-curve slopes for slender bodies of revolu- 
tion. For example, from Appendix C there is obtained 


d% 

<1°' a=o 
where 


= 8 ( ^ 


1 7+1 

.cJ 


( 39 ) 


*3-Cn 

da 


TCN 

MsV-i 

Mk. 


1 + 2|3 1 + ^ 


1+7(MoSn)‘ 


-ig; 


(Mo&n)‘ 


— 

7+1 \°s. 


1 + (MoSu)^ 


Hence, the pertinent flow parameters necessary to determine the initial 
slope of the normal-force-coefficient curve by means of the above expres- 
sions may be obtained from table I for the case of a/Sjj = 0. 


EXPERIMENT 


In order to obtain a check on the predictions of the preceding 
theoretical analysis, the pressures acting on the surfaces of lifting 
bodies of revolution corresponding to values of the hypersonic similarity 
parameter K from 0.60 to 1.68 at Mach numbers from 3-00 to 5*05 were 
determined experimentally. The bodies were tested at angles of attack 
up to 15*^- A brief description of these tests follows. 


Test Apparatus 


Tests were conducted in the Ames 10- by l4-inch supersonic wind 
tunnel. A detailed description of the wind tunnel and auxiliary equip- 
ment may be found in reference 10. The pressures acting on the model sur- 
faces were measured with a mercury U-tube manometer or by means of McCleod 
gages when the pressures were low enou^ to be recorded on the latter. 

Pressure-distribution models were mounted on a 0^ model support and 
on 10^, and 15^ bent supports. The test models were two tangent ogives 
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having fineness ratios 3 and 5 and two cones having the same vertex angles 
as the ogives. The dimensions of these models and location of the pres- 
sure orifices are shown in figxjre 3* 


Tests and Procedure 


Pressures on the model surfaces were measured at 0*^^ 10*^^ and 

15° angles of attack at test Mach numbers of 3-OOj ^-25, and 5 •05* The 
Reynolds numbers (based on maximum diameter of the ogives) were 1.09 
million at Mach numbers 3-00 and 4.25, and 0.52 million at Mach number 
5 . 05 . 

The pressures around the cone surface (0° to 3^0°) at meridian sta- 
tions 45° apart were recorded simultaneously at each Mach number and 
angle of attack. In the case of the two ogival models^ the pressures 
were recorded at meridian stations 9C apart^ Each model was then rotated 
45° about its longitudinal axis (except at 0 angle of attack) and the 
process repeated. 

Schlieren photographs of the bow shock waves in three meridian planes 
were also obtained. 


Accuracy of Test Results 


The variation in Mach mamber from the nominal value did not exceed 
±0.02 in the region of the test section where the models were located. 

The precision of the computed pressiire coefficients was affected by 
inaccuracies in the pressure measurements, as well as uncertainties in the 
stream angle and the free-stream dynamic pressure. The resulting errors 
in the pressure coefficients were generally less than ±0.005 throughout 
the Mach number range for all angles of attack. The meridian angles at 
which the pressure coefficients are plotted are considered accurate to 
within ±1°. 


COMPARISON OF THEORY WITH EXPERIMENT AND DISCUSSION OF RESULTS 

Flow at the Vertex 


It will be recalled that one of the fundamental assumptions in the 
(3^gY0Popment of the conical flow theory was that the conical shock remains 
circular when the cone is inclined. It is appropriate to examine the 
validity of this assumption before proceeding with a comparison of the 
theoretical and experimental surface pressures. To this end, schlieren 
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evidence on the shapes of the conical shocks for the twQ test cones at 
Mach number 5*05 is presented in figures 4 and 5 fo^r angles of attack of 
0^ , 5^^ 10^^ and 15^. The data shown in the figures were obtained by 
measuring the angle of the shock in the schlieren photographs at various 
meridian stations. Theory is compared with experiment in a cross-sectional 
plane at an arbitrary distance downstream of the vertex. It is observed 
in figure 4 that the conical shock attached to the slender cone remains 
nearly circular for angles of attack up to 10^. At a = 15^, the angle 
of attack is greater than the cone half-angle and^ as mi^t be expected^ 
the shock is no longer circular. Hovever, in the case of the blunt cone 
(fig. 5 )^ "the conical shock remains essentially circular for all angles 
of attack up to and including 15*^. It would appear, then, that at least 
for moderate cone angles, so long as the angle of attack is less than the 
semivertex angle of the cone, the assumption of a circular conical shock 
made in the analysis is justified. 

The second basic assumption employed in the development of the conical 
flow theory of this paper is that co - 6 < < 1. It is apparent that this 
condition is best satisfied for blunt cones and for hi^ Mach numbers. 

The accuracy of the theory might be expected, therefore, to improve both 
with increasing cone angle and increasing Mach number. The predictions 
of theory and the results of the pres sure -distribution tests for the two 
test cones (5c = 11.42^ and 6c = l8.92^) are shown in figures 6, 7, and 8. 
The data are plotted in the form of surface pressure coefficient as 
a function of the meridian angle cp.® It is observed in these figures 
that the predictions of theory, when applicable, are in good agreement 
with experiment for the Mach numbers and angles of attack presented. It 
is evident also that at the highest angle of attack (a = 15 ^^) the theory 
is less reliable on the leeward side of the body. Althou^ this result 
is due in part to the limitations of the theory, it is also clear that 
the viscous effects of the flow are influencing the pressures to a greater 
extent over the leeward portion of the body. It can also be deduced from 
these figures that agreement between theory and experiment is better for 
the blunter cone, particularly at the hi^er angles of attack. It is 
indicated, therefore, that the predictions of the conical flow theory of 
this paper will yield more reliable results when the parameter a/s^ < 1. 

In the lower range of Mq and 5c (fig. 6) where the present conical flow 
theory is not applicable (see fig. 2), Stoners second-order solution 
(ref. 11 ) applied in the manner described in reference 12 yields results 
which are in good agreement with experiment.^ It is observed in figure 7 

^It will be noted in these and subsequent figures that the data are 
often plotted at meridian stations slightly different from 0^^ 45^, 90 *^, 
etc. This resulted from inaccuracies in rotational positioning of the 
model . 

"^Due to the limited range of the results presented in the tables of 
reference comparison can be made only for the slender cone and then only 
for Mach numbers 3-00 and 4.25 without resorting to extrapolation. Although 
the agreement between Stoners results and experiment appears to be better 
at a = 10*^ than at a = 5^^ this result must be attributed to the manner 
in which the flow parameters presented in reference 5 were interpolated. 
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that both methods yield comparable results over the angle-of-attack range 
at Mach number 4.25« Hence ^ it appears that the two theories tend to 
overlap in regards to their usefulness over the range of Mq and Sc- 

It should be mentioned here that the treatment of inclined cone flow 
as presented by Ferri in references 7 ^•iid I3 was applied to the cases 
under discussion in the present paper. Ferri* s method did not yield 
results as good as either Stone's second-order theory (where applicable) 
or the conical flow theory of the present paper. In fact, inconsistent 
results were obtained when the method was applied according to references 7 
and 13. This discrepancy may be traced to what appears to be an inconsist- 
ency between equation (55) in reference 7 and equation (55) in reference I3. 

From the preceding comparisons of the experimentally determined sur- 
face pressures and shock-wave shapes with the predictions of the conical 
flow solutions of this paper, it is indicated that the latter solutions 
may be employed to predict the flow properties about a lifting cone at^ 
hif^ supersonic airspeeds with good accuracy when the angle of attack is 
less than the cone half-angle (i.e., when cc/Sq < l) . It is therefore^ 
suggested that these solutions may be particularly useful for determining 
conical flow fields about lifting cones over the range of Mq and not 
treated in the M.I.T. tables (ref. 5). 


Flow Downstream of the Vertex 


It remains now to determine the accuracy with which the solutions 
for flow about lifting cones in combination with the isentropic expansion 
equations predict the flow about bodies of revolution other than cones. 

The pressure distributions on the surfaces of two ogives (having fineness 
ratios 3 and 5) at Mach numbers 3-00, 4.25, and 5.O5 and at angles of 
attack of 0°, 5°, 10°, and 15° were calculated using the methods of this 
paper. These distributions are presented in figures 9, 10, and 11 for 
values of the hypersonic similarity parameter, K, varying from 0.60 to 
1.68. Also shown are the results of the pressure-distribution tests for 
the two ogival models . 

Comparing first the predictions of theory with experiment for the 
case of zero lift (fig. 9), we observe that the accuracy of the shock- 
expansion method generally improves with increasing K. This trend is, 
of course, the same as was observed in reference 1 for comparisons of the 
predictions of the shock-expansion method with those of the method of 
characteristics. The results of a characteristics solution for the 
Z/d = 3 ogive at Mq = 3. 00 (from ref. l4) are also shown for comparative 
purposes. Characteristics solutions are not available for the other cases; 
however, the results of Rossow (which were obtained by correlating the 
pressures obtained by characteristics solutions according to the hypersonic 
similarity law and presented in ref. 15) are shown. As mi^t be expected, 
Rossow* s results are generally in good agreement with experiment althou^ 
there is a sli^t imderestimation of the pressures near the base of the 
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body at Mq = 5 •05* This is attributed to the viscous effects of the flow 
which are probably influencing the pressures at this Mach number. In any 
event, it is evident that the predictions of Rossow and the shock-expansion 
method are in good agreement at the highest value of K (K = 1.68). 

It is of interest now to determine the reliability of the predictions 
of the shock -expans ion method for lifting bodies. As shown in figure 10, 
the theory yields good agreement with experiment for the fineness-ratio-5 
ogive on the windward side of the body except at Mq = 3.00 (K = 0.60).® 
Some disagreement is evident, however, on the leeward side of the body 
at all Mach numbers. In the case of the fineness -ratio-3 ogive (fig. ll) , 
agreement is generally better over the entire body at each angle of attack, 
particularly at the higher values of K (K > l) . It will be recalled from 
figure 9 that at a = 0° the longitudinal pressure distributions on both 
ogives indicated that the accuracy of the shock-expansion method increased 
as K increased. It is indicated in figuies 10 and 11 that, as would be 
expected, this trend carries over to the case of lifting bodies. 

It appears in figures 10 and 11 that the most important factor 
influencing the accuracy of the method is the reliability of the conical 
flow theory at the vertex, since the inaccuracies at the vertex appear 
to be reflected strongly in the pressures downstream of the vertex. The 
question naturally arises, then, how good are the predictions of the 
shock-expansion method when experimentally determined initial conditions 
at the vertex are employed? To answer this question, the pressure coef- 
ficients on the surfaces of the two ogives under discussion were deter- 
mined in the following manner. Initial conditions at the vertex were 
determined from the measured static pressures around a cone (corresponding 
to the vertex angle of the body) in conjunction with the measured shock- 
wave angle (in the plane cp = 0) obtained from schlieren photographs of 
the conical flow field. The pressure coefficients downstream of the ver- 
tex were then calculated as before. The results of these calculations 
for Mach numbers 3.00, 4.25, and 5-05 are compared with experiment in 
figures 12 and 13 for a = 15°. Results for a = 15° are presented because 
at this angle of attack the applicability of the conical flow solutions 
is most marginal. It is observed from figure 12(a) that the theory yields 
results which indicate an underexpansion of the flow on the sides of the 
body (cp = 45° and cp = 90°). This result is not surprising since, at this 
low value of K (K = 0.60) , it would be expected that the true streamlines 
would deviate considerably from a meridian line. In other words, flow 
disturbances in planes tangent to the body at the surface are no longer 
small compared to those in axial planes. Hence, the flow along a true 
streamline travels through a greater resultant angle than that represented 
by a meridian line. It can be seen from figures 12(b) and 12(c) that as 
the Mach number and, hence, K, is increased, better agreement is obtained. 

°It should be noted in figure 10(a) that Stone’s second -order solution 
is employed at the vertex since the conical flow theory of the present paper 
is not applicable for these conditions (i.e., Mq = 3-00 and 8c = 11.42°; 
see fig. 2) . 


26 


MCA TN 33^9 


This is attributed to the fact that the streamlines of the flow deviate 
less from meridian lines as K is increased. The same general trend may- 
be noted in figure I3 for the fineness-ratio-3 ogive. However, in this 
latter case, over-all agreement appears to he better. In fact, good 
results are obtained for values of K > 1 except on the extreme leeward 
side of the body where it is probable that viscous effects are influencing 
the pressures. There may be some separation of flow over this portion of 
the body although no evidence of this could be determined from the schlie- 
ren photographs. In the case of the fineness-ratio-5 ogive, schlieren 
evidence indicated flow separation on the leeward side of the body for all 
Mach numbers at a = 15 It may be deduced from these figures that the 
application of the shock-expansion method will yield better results when 
the initial conditions at the vertex are determined from cone tests 
rather than from presently available cone theory. 

There now remains the determination of the accuracy of the predic- 
tions of the generalized shock -expans ion method for the flow field (other 
than the surface) about a lifting body of revolution. To this end, flow 
in the plane of symmetry (cp = 0° and cp = l80°) was calculated (after the 
manner discussed in Appendix B) for each ogive traveling at a Mach number 
5.05 and at an angle of attack of 10°. ) Flow in a side meridian plane 
(cp = 90°) was also calculated for the fineness -ratio-3 ogive. The result- 
ing shock-wave shapes are compared with the actual shapes (obtained from 
schlieren photographs) in figiore l4. The theoretically determined conical 
shocks are also shown for contrast. In the case of the fineness-ratio-3 
ogive (K = 1.68 and a/5N = 0.53 ) > theory and experiment are observed to 
be in excellent agreement in the plane of symmetry. The same observations 
may be made for the side meridian plane. In this latter connection, it 
is of interest to point out that essentially the same result would have 
been obtained if the shock were assigned circular in cross-sectional planes 
and its location determined from the calculations in the plane of symmetry. 
In view of the agreement between theory and experiment, it is indicated 
that when K is greater than 1 and a/b^ is about 1/2 or less, the shock 
is essentially circular in cross-sectional view of the flow field about 
a pointed body of revolution. In the case of the fineness-ratio-5 ogive, 
the poor agreement on the leeward side of the body might be expected since 
not only is K marginal for the application of the theory but, more 
important, a/5jj is relatively large (0/5^ = 0.88). It should he pointed 
out that if experimentally determined initial conditions are employed, 
good agreement with experiment downstream of the vertex is obtained. 

Although the predictions of the generalized shock-expansion method 
have been checked only at the inner and outer boundaries of the flow 
field, it is expected that equally good results would be obtained at 
intermediate points in the flow field. This conclusion is based on the 
fact that the bow shock waves were obtained as a result of the calcula- 
tions of these intermediate points. 
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Hypersonic Slender-Body Theory 


The predictions of the hypersonic slender-body theory for a = 0° 
and a = 5 ° are compared with experiment in figures 15 and I 6 , It appears 
from a comparison of figures 9 and I 5 that the slender-body theory will 
yield more accurate drag coefficients than the more general theory at 
a = 0°, particularly at the lower values of K. This result, although 
doubtlessly fortuitous, is the same as that found in reference 1. In the 
case of lifting bodies, it appears that the slender-body theory yields 
results which are somewhat less satisfactory at all values of K. However, 
the theory displays sxifficient accuracy for many engineering purposes 
even at K = 1. This is particularly evident for the more slender of the 
two bodies as indicated in figure l 6 (b) . In view of its simplicity, the 
slender-body method should prove useful and its application is further 
facilitated by the presentation in this paper of tabulated values of the 
pertinent flow parameters for selected values of MqSu and a/Sjf (see 
table I), 


Normal-Force Coefficients 


It is appropriate now to consider briefly the forces experienced by 
the ogives. To this end, normal-force coefficients were obtained by 
integrating the theoretical pressure distributions for the two ogives at 
a Mach number of 5* 05- The results of these calculations are compared 
with those obtained from integrated experimental pressure distributions 
in figure I 7 for values of K of 1.01 and 1.68. It is observed that 
theory yields values of Cjj which are, in general, higher than those 
obtained by experiment. However, agreement improves with increasing K. 
Equation ( 36 ), as well as the hypersonic slender-body solution (eq. (39)), 
appears to yield satisfactory initial normal -force -curve slopes at values 
of K as low as 1. Axial forces have also been obtained for these ogives 
and the shock-expansion method is found to apply with essentially the same 
accuracy. 


CONCLUDING REMARKS 


The flow about a lifting body of revolution at high supersonic air- 
speeds was investigated analytically. With the assumptions of conical 
flow at the vertex, hi^ supersonic Mach numbers, and low angles of attack, 
simple approximate solutions were obtained which yield the Mach number 
and pressure distributions on the surface of the body. Surface stream- 
lines were approximated by meridian lines and the flow field in meridian 
planes was calculated by means of the generalized shock-expansion method. 

In the special case of slender bodies, simple explicit expressions were 
obtained for the Mach number and pressure distributions on the surface. 
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Surface pressures and shock-wave shapes were obtained experimentally 
at Mach numbers from 3-00 to 5-05 for two ogives having fineness ratios 3 
and 5 for two cones having the same vertex angles as the ogives. The 
predictions of the methods of this paper for the surface pressures were 
found to be in good agreement with experiment at values of K about 1 ^ 
or greater^ when a/b]^ (the ratio of angle of attack to semivertex angle) 
was less than about 1 / 2 . For increasing values of this parameter, agree- 
ment deteriorates but may still be considered fair for values of 
up to about 1 . The generalized shock-expansion method yielded very good 
agreement with experiment for the shape of the bow shock at K = 1.68 and 
= 0 . 53 - Tt was further indicated that the bow shock remains essen- 
tially circular in cross section for angles of attack up to approximately 
one half the semivertex angle when K is greater than 1 . 
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APPENDIX A 

DETERMINATION OF FLOW CONDITIONS AT THE VERTEX IN THE PLANE OF SYMMETRY 


Due to the symmetry of the flow in the plane of symmetry. 


^ = 0 


Sep 


= 0 w = 0 


Now from the flow geometry 


u = V cos(o) - 6) 

V = -V sin(co - 6) 

and from equation (la) 

Su 

0(jl) 

The flow is therefore irrotational and the following relation holds 


V “ 


- tan(co ^ 5 )dS 


(Al) 


Substituting the above expressions in equation ( 5 ) results in 
2 


7 ■ 1 


cot(o) - 5)( 1 -I- - cot CO t tan(co “ 


csc(g) - 5) Sw 
V sin 0 ) ^cp 


- tan(w ~ 1“ = (A2) 


Equation (A2) is not amenable to algebraic solution. However, since it 
differs only in the term containing bw/brp from the equivalent equation 
defining the axisymmetric flow field, a solution analogous to that employed 
in reference 1 is suggested. Consistent, then, with the restrictions 
imposed on the flow field in this analysis, namely, (w - S) < < 1 radian, 
equation (A2) reduces to 
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1 - (o) - S)cot (J + [1 - M^(o) - 5)2] i ^ = 0 (A3) 

c3u) V sin 0) dcp 


vhere M and^ hence^ the magnitude of V are considered constant (see 
ref. 1). Now from equation (15) 


3w sin (jOo 

— = ±V e ^ 

dcp ° sin 0) 


(1 


- n 


cot lOg) 


(a4) 


where dw/3cp is positive at cp = 0 and negative at cp = «. Hence, near 
the surface of the cone, equation (A3) reduces to the linear equation 




-(1 + a) 


or 


(jO - S = (w - 6c) (2 + a) 


(A5) 


where (since q < < l) 


a 



sin u)g 
sin^&c 


Combining this relation with equations (A3) and (A5) results in 


^5 (u) - &c)(2 + a) cot w - 1 o sin^6c 

Su ■ 1 _ m2(u _ 6c)^(2 + a)^ sin2w[l - m2(w - &c)^(2 + 0 )^] 

(a6) 


which can be integrated to yield (substituting in the boundary conditions) 


5 - 6c = a sin^5c(cot u) - cot 6c)^l - ^ sin^&c|csc® 


5n - 


M(2 + a) 


+ CSC 


5c + 


M(2 + a)J 


} 


(2 + a) taxi 6csec^6c 
(2 + a)^M^tan^6c - 1 


ln[l + (w - 6c) cot 6c] + 


sin‘ 


6c - 


M(2 + a)_ 


•|tan 6c + M[1 + a + (2 + a)M taxi 6cjJ - a sin26cM[l - (2 + a)M tan 6cJ 


2 sin" 


M(2 + a)_ 


(2 + c)M [1 - (2 + c)M tan 6c] 


ln[l + (2 + ct)M(w - 6c) ] + 


sin^ 

L M(2 + a)J 

■jtan 6c 

+ M[(2 + cr)M tan 6c - (l + 0 )]^ + oM sin^6c[l + (2 + c)M tan 6c] 

2 sin^ 

5c + ““7 ^ 

M(2 + a)_ 

(2 + 0)M^[1 + (2 + a)M tan 6c] 


ln[l - (2 + cr)M(w - 6c)] 


(AT) 


If M in equation (AT) is taken as the Mach number just downstream of the shock, then flow 
conditions at the shock (i.e., Ms, tOs^ sxid 6s) can be obtained by the simultaneous solution of equa- 
tions (6)^ (T) ^ (at) ^ and the oblique shock-wave relations 
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[ 27 Mo^sin^((jJs ± a) - (7 “ l) ] [ (7 ~ l)Mo^sin^((jg + a) +2] 


(A8) 


and 


Bg ± ot' = tan"^ 



(A9) 


To apply these equations, a value for e is chosen slightly less than 
a. Then q is fixed by equation (6) and the relationship between the 
shock-wave angles is obtained by setting cp = jt in equation (7), namely. 


Now (wg)^_Q is determined by solving equations (A7) , (A8) , and (A9) 
simultaneously where a > 0 in equation (A7) and the positive sign is 
used in equations (A8) and (A9) , Then (ti)s)q)=jt is determined in the 
same manner as (ws)m=o where now a < 0 in equation (A7) and the negative 
sign is used in equations (A8) and (A9) . If the resulting value of e 
(calculated from eqs. (6) and (AlO)) differs from that originally chosen, 
the procedure is repeated using the calculated value of e, and so forth. 
Althou^ the foregoing procedure is somewhat tedious, the number of itera- 
tions can be reduced to two or three in most cases by carefully choosing 
e and (u)s)q)=o- this connection, it has been found useful to choose a 

value for (103)35=0 which is less than the corresponding shock-wave angle 
of the nonlifting body by q , the latter angle being approximately 10 to 
15 percent of a at the hi^er Mach numbers (Mq >4) and I5 to 30 percent 
of a at the lower Mach numbers (Mq < 4) . It should be remembered that 
K is always approximately 1.4 or greater in this analysis. 

It is clear that equation (A7) should give a better representation 
of the flow field in the cp = 0 plane than in the cp = rt plane since 
(co - 5) is always less on the hi^ pressure side of the body. In fact, 
there are cases when the combination of Mq, and a is such that 

equation (A7) no longer applies on the leeward side of the body. For 
example, if the body is slender such that the angle of attack approaches 
the half -cone angle, |a| can be 1 or greater. Since a < 0 in this half- 
plane (cp = it) , 8s/8(j > 0 for these conditions and equation (A7) will no 
longer represent 5 as a monotonically decreasing function of u). Hence, 
it is possible that no simultaneous solution of equations (A7) , (A8) , and 
(A9) will exist in the cp = it plane. It is necessary, then, to obtain 


(ws)(p=jt " (‘*^s)cp=o " 2q 


(AlO) 


another expression relating 8 and w in this plane. The development of such a relationship by- 
imposing the restriction that 8 < < 1 radian but u remains arbitrary will now be considered.^ 


With the restriction that 8 < < 1 radian, equation (A2) may be reduced to the form 

8 cot u + ^ [1 - M^sin^cj(l - 6 cot w) ^ = 0 

oo) V sin to ocp 

where M is again considered constant. Near the surface of the cone, the above expression 
to 

^ + 8 cot u) + a = 0 
dto 


which has the solution 


8 = a cot to + k CSC to 


Combining equations (All) and (A12) and integrating the resulting expression yields (to the 
of accuracy of this analysis) 


_ 8 _ 

Sc 


1 + 


1 + 


2M^(a+ l) (sin^Sc - sin^to) 

1 + *i(l+ oM^sin^8c)^+ UM^sin^Sc - sin^Sc 

2M^(g+ l)(sin^8c - sin^to) 

1 "^/(l+ oM^sin^8c)^ + i(M^sin^8c “ aM^sin^Sc 


M sin 5c 


2n/(i + aM^ sin^8c)^ + 4M^sin^8c 


[2 - sin2to(l + a csc^Sc)]^ ' 

] "l^gin 8 (2 + o)M^sin^8c - (l + a)M^sin^to] -M^sin^8c^ 


a(a + sin^? 


\2a -M^sin^Sj. 


y[2 - (a + sin^Sc) 


(All) 

reduces 

(A12) 

order 


a^) 

(A13) 


^The method of development is similar to that presented in reference 1. 


Lx) 

Lx) 


NACA TN 331^9 


34 


MCA TW 3349 


where M is again chosen as the Mach number just downstream of the shock. 
Hence, conditions at the shock can be determined by the iteration pro- 
cedure mentioned earlier except that now equation (AI 3 ) is used in place 
of equation (A?) . Althou^ the former equation may be employed to deter- 
mine the shock angle at cp = 0 as well as at cp = ir for slender bodies, 
it is suggested that equation (AT) be used to determine flow conditions 
on the windward side of the body and equation (AI 3 ) for flow conditions 
on the leeward side of the body only when equation (AT) cannot be solved 
simultaneously with equations (a 8 ) and (A 9 ) . If the shock angles in the 
plane of symmetry, are known, the variation of 5 with cj is known from 
equation (AT) . For the case where this equation does not apply in the 
plane cp = ir, this variation is given by equation (AI 3 ) . 


The determination of the small variations in the local velocity in 
the plane of symmetry is identical to that presented in reference 1 for 
the case of (u) - 6) < < 1 radian. Hence, the expression for the velocity 
may be written 


Vs 


5 c- 5 e 


cos(ws - 5s) 




■(5c - 5)( 


U)c 


5s) 


-> 


cos 


5q “ ^s 


(Al 4 ) 


If the velocity is known, the Mach number may, of course, be deter- 
mined from the relation 


M^= 



(AI 5 ) 


and the pressure coefficient may be obtained with the aid of the expres- 
sion 


P 2 fVs V 

7 Mq^ \Po Ps 



(AI 6 ) 


where 

Pg 27Mo^sin^(cJs ± o.) - (7 - 1 ) 

% " y 1 


and 
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APPENDIX B 

DETEIMINATION OF FLOW IN A MERIDIAN PLANE 


Flow about a lifting body of revolution may be calculated in any 
meridian plane by the generalized shock-expansion method in much the same 
manner as the procedure employed in reference 4 for flov about airfoils. 
However^ the application of the method is somewhat more complicated for 
the case of a body of revolution since now the influence of the conical 
flow in the region of the vertex must be considered. 

The determination of axially symmetric flow in the region of the 
vertex of a body of revolution (K > l) was described in reference l6. 
Expressions were developed which yield the shock-wave curvature as well 
as flow conditions along a line a short distance downstream of the vertex. 
An analysis entirely analogous to that in reference l6 was carried throu^ 
for the lifting body and it was found that more general expressions can be 
obtained which take into account the effects of angle of attack. Thus, 
it can be shown that the expression for the ratio of the shock-wave cur- 
vature to body curvature near the vertex is (consistent with the assump- 
tions of the present paper) 


Kg 27M0 sin^[o)s +a(cos cp + T] cot cog sin^cp) ] - (7 - 1) 
% 


— g ; — Tl -tan(u)s -Su)cot 

Mo sin'^Piysin 2[us +o.(cos cp + T) cot WgSin^cp) J <• -■ 


where 


(Bl) 


‘*>8 = (‘^s)(p=o + ■ cos cp) 


The function ^ is defined by the expression 


(l+ 1 


V V sin 5jj dCp/ 

L vv J 


and is evaluated at the surface outside the vortical layer by means of 
the previously developed conical-flow expressions. Similarly, expressions 
for flow conditions along a line normal to the axis of the body a short 
distance downstream of the vertex may be obtained. For example, the 
variation of flow inclination, 5, along this line is given by 


5=5a- 


sin Sacos &a 

ya 


(y -ya) 


5b - 5a + 


\(f sin 5 acos 5a 
7A 


(yb-Ya) 


y - Ya 

yB-y. 


(B2) 
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vhere A and B are points on the line corresponding to the surface and 
the shocks respectively. The relations for the static pressure and total 


respectively. It should he noticed that expressions (Bl) and (B 2 ) reduce 
identically to those given in reference 16 for axially symmetric flow. 

Knowing the flow .conditions in any meridian plane in the region of 
the vertex, it is now a relatively simple matter to construct the entire 
flow field downstream of the vertex. To illustrate, consider the sketch 
(flow in a meridian plane cp = constant): 


pressure remain unchanged and may he written (in the notation of the 
present paper) 



2 tpK]^cos(P -f 5 ) 
sin p sin 23 



2 


(B3) 


and 






y - ys 

- y 



2 


Ha - Hb - 


(b 4 ) 


y 



X 
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With the ohlique-shock-wave, conical-flow, and expansion equations, all 
fluid properties at points M, A’, A, C, and so forth, on the body surface 
may be calculated in the manner described previously in the present paper. 
Hence, flow conditions along the line AB may be determined from expres- 
sions (Bl) through (b 4 ). It will be recalled that the basic condition 
employed In constructing flow fields about airfoils by the generalized 
shock-expansion method is that the pressure is constant along Mach lines 
emanating from the surface. In the case of flow about bodies of revolu- 
tion, this condition must be relaxed to account for the variation in pres- 
sure due to the influence of the conical flow in the region of the vertex. 
This may be accomplished in the following manner. The Mach line A'B is 
constructed from the known conditions in the region MAB shown in the 
sketch. The net pressxire change along this Mach line (i.e., pg - p^i ) 
is thus determined. This pressure difference is then assumed to represent 
the net pressure change between the body surface and the shock along each 
Mach line emanating from the surface downstream of the vertex. The flow 
field is then constructed using this criterion in conjunction with the 
isentropic expansion relations for flow along streamlines. Once the 
shapes of the streamlines are calculated, the fluid properties along these 
lines are, of course, determined in the same manner as those along the 
surface. 
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APPENDIX C 


FLOW AT THE SURFACE OF A SLENDER BODY TRAVELING AT HIGH 
SUPERSONIC AIRSPEEDS AND AT SMALL ANGLES OF ATTACK 


If a slender body (i.e,, a body on the surface of which the slopes 
are everywhere small compared to l) is traveling at free-stream Mach 
numbers very large compared to 1 (again, of course, K > l) and at angles 
of attack very small compared to 1, the local Mach numbers will likewise 
be large compared to 1. It follows, then, that the inclination of the 
nose shock wave will be small and, consequently, that w will always be 
small. In this case, the relation between 5 and u) at the vertex (in 
the plane cp = O) is relatively simple and may be obtained by integrating 
the expression (see eqs. (Ak) and (All)) 


8S _ 5 €^s 

bji ~ ~ ^ ~ 0)2 

which yields 



(Cl) 


Combining this expression with equation (Al) , the relation 




(C2) 


defining the velocities in the plane cp = 0 is easily obtained. Hence, 
the surface Mach number in this plane at the vertex, may (to the 
order of accuracy of this analysis) be related to Ms by combining 
equation (AI 5 ) with equation *(C2) to yield 


m/ 



(Ms5n) 





e ^ 



(C3) 


Now the oblique shock-wave relations for flow of the type under 
consideration reduce to (at cp = O) 
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Ms^ = 


(7 + l)^Mo^(wg + a)^ 


[27Mo2(tos + a) - (7 - l)][(7 - 1 )Mo^(u;s + a)^ + 2] 


(C4) 


Mo^(cjs + a)^ - ^ Mo^(u)s + a) (S^ + a) -1 = 0 


( 05 ) 


and 


Ps 

Po 


27 

7 + 1 


Mq ( 


+ a) = 


7-1 
7 + 1 


Combining equations (6), (Cl), and (C5) results in 
MQ(ws+a) 


(C6) 


f /(Moa+^MoTi ) +3 


1+^^i^o^ ( 6W+a) ( 5N-n ) 



There remains now the determination of ti, which defines the position 
of the conical shocks in order to determine the shock-wave angle in the 
plane cp = 0. To this end, the assumption of a circular conical shock is 
again employed, but now it is deemed sufficiently accurate for the purposes 
of this analysis to assume a linear variation of t) with a; namely 


^ = 


da 


a 

a=o 
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Now from equation (AlO) there results 


dq _ 1 / du)s 
dS; " 2 \ da 

Consider for the moment, flow in the plane cp = jt. The angle of attack 
a and, therefore, e are both negative in this plane. Hence, equa- 
tion (C5) may he written 


Mo^(o.'s - a)^ - Mo^(tOs - a) (Sg - a) -1=0 

and the conical-flow expression (see eq. (Cl)) becomes (at co = ug) 


cp=rt 


dojg 

da 


cp=o 


(C9) 


SgWg = + 6 (los " hN) 


Differentiating these expressions with respect to a and combining the 
resulting expressions yields in the limit as a — ^ 0 


2u)g( 


'dcog 
, Sa 



7 + 1 
2 


^ (u)g - 5 n) 
da 


(us + Ss) 


= 0 


(CIO) 


Preceding in the same manner for flow in the plane 9=0 (see eqs. (Cl) 
and (C5)) there is obtained 




7 + 1 


ha 


((jjg - 5^) + (o)s + 5s) 


= 0 


(Cll) 


Combining equations (6) , (C9) , (CIO) , and (Cll) results in 


dq 


da 


a=o 


1 - 


1 + 


y + b 

^ \ ‘*>s y 


y + b 4jg - &N 

4 V wg , 


At a = 0 (see ref. l) 
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4l 


5N 


and 


Mo“s ^ + 2 (MqSn) 

Hence, 


MqSjj 


1 + 


7 + 1 


(Mo5n)‘ 


fC12^ 


11 


1 


7 + 1 ^ 

6n V 



1 + 


7 + 1 

~T~ 




a 


(CI3) 


where ^s/^N is a function of Mo5]\f only and, of course, is given by 
equations (C12) . The shock-wave angle at cp = 0 can now be determined for 
given values of Mq5j^ and Moa by means of equations (C7) and (CI 3 ). 

Hence, the shock angle around the entire conical shock front may easily 
be determined with the aid of the expression 


tJs = (“s)q,=o “ cos cp) (Cl4) 

Surface conditions around the body at the vertex may now be deter- 
mined after the manner described in the more general analysis of flow 
about cones. For example, consistent with the assumptions basic to the 
present analysis, equations ( 17 ) and (I 8 ) reduce to 


'W'e __ (£s 

(Vn)cp=o ” 


sin cp 


(CI 5 ) 


and 


^e 

(%)cp=o 


= 1 + e(jJs(l 


cos cp) 


(C16) 


respectively. Hence, the surface velocity external to the vortical layer 
at the vertex may be written 


42 


MCA TN 3349 


2 


~ Ve 
J^n)cp=o. 


= 1 + 2eus(l - cos 


cp) 



(CI 7 ) 


and the corresponding Mach number. Mg, may be related to (M|^)(p=:Q by 
combining equation (I 9 ) with this expression to yield 


( Mk)(P=o 

Me 


= 1 - 


7 - 1 




2 (1 


COS cp) 




(CI 8 ) 


Consistent with the assimiptions basic to the present analysis, the Mach 
mmiber directly on the s\irface at the vertex (i.e., inside the vortical 
layer) may be obtained from the relation (see eq. (20)) 


where 

Ec-E 

7Cv 


e 


e 



(CI 9 ) 


— — 1 

27^ ^)cp=0 "■ (7 " 1) 

1 

7 

Mq^(ws + O' cos cp)^ \{7 - 1 )Mq^(ws + o.)^=o + 

27Mo^(ws + <x cos cp) - (7 - 1) 


Mo^(ws + a)^ f(7 - l)Mo^(oJs + acos cp)^+ ^ 

k cp — 0 ^ t 


Equations (C3) , (Cl8), and (CI 9 ) may now be combined to yield 


(C20) 


7 - 1 

(Ms)J=o^l+— (MgSw)^ 




1 + 


- (|f) 




Ec~Ee 
> 7Cv 






COc 


2— (l-coscp)+v5j^ 


5i) -i-'"'!’]} 


(C21) 


The expressions just derived provide the Mach number on the surface 
at the vertex. If M]^ is known^ the Mach number anyvhere on the surface 
of the body may be obtained by means of the expression (see ref. l) 



(C22) 
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Now the pressure coefficient is given by the expression 


P = 




Ps 


(Ps)q)=o 


- 1 


(C23) 


The pressure rise across the shock is given by equation (c6) and the ratio 
of the pressure anyvhere on the surface to the pressure at the shock at 
cp = 0 can he expressed (to the order of accuracy of this analysis) in 
the form 


R 


(Ps) 


cp=o 


■(Ms),p=o' 


M 


7-1 


This expression may be combined with equation (C22), and equation (C23) 
cau be written 


P 


2 

7Mq^ 




M, 


•N 



(MsSn)cp=o 



^ 27 



1 


(C24) 


yielding the pressure coefficient at any point on the surface of the body. 

The initial normal -force -curve slope for slender bodies of revolution 
may easily be determined in the following manner. To the order of accuracy 
of this analysis^ equation (36) may be reduced to the form 


dC^ 


da 


0=0 





dCN 

W IT 


TCN 



Pn 

PN P 



(C25) 


Consistent with the assimptions basic to the present analysis^ the fol- 
lowing relations may be obtained; namely. 
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and 


JL 

pn 



These expressions may be combined with equation (C22) and (C25) to yield 


d% 

da 


a=o 



r\ 

) 

TCN J 
0 




111 


7^7 (C26) 


There remains now the determination of the initial normal-force-curve 
slope for a cone tangent to the body at the vertex. This slope may be 
expressed as 


^ ^ T” A 

da tcN da 

NoW the ratio of the static pressure anywhere on the surface of a lifting 
cone (a < < 1 ) to the free-stream static pressure may he expressed in the 
form 


j cos cp dcp 


(C2T) 


% 

Po 



(C28) 


where Ve is the local velocity externally adjacent to the vortical layer 
and 


H 

■ (7 + l)Mo^(us + a- cos op)^ “ 

7 -x 

2 2 n 

27M0 (wg +a cos cp) - (7- 1) 

Ho 

- 1)Mo^(us +“• cp)^ + 2_ 


y + 1 


(C29) 


Differentiating equation (C28) with respect to a and retaining only 
terms which sire functions of cp, we obtain in the limit as a — > 0 

f \ 



% y 

Po 





(Mo6n> 


(MqWs) 


2Al 


1+^ (MqSn) 


1 + 7(MoSn)' 


(C30) 


A 

da 
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From equations (6) and (CI3) there may be obtained 


de 

da 



Combining this expression with equation (C30) and noting that (see ref. l) 



[1 + 7(Mo5n)^] 


MH-y 


27 

7-1 


equation (C27) may finally be written 


dCjNj 


da 


TCN 


27 


yMsV-" 

w 

1+^(Mo5n)^ 

{ 

1+^(MoSn)^ 

1+7(Mo5k)^ 


(MqSn) 

7+5 5k 
7 + 1 5s " _ 

1+ 2 " (Mo&n) 



(C31) 


where Mg, and 5g may be determined from reference 1 or by setting 

a = 0 in the pertinent expressions previously derived in the present 
paper. 
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TABLE I.- TABLE OF FUNCTIONS FOR HYPERSONIC SLENDER BODY METHOD 



\a/SN 

0 

0.20 

0 . 4 o 

0.60 

0.80 

1.00 

0. 60 
.80 

1.00 

1.20 

1. Uo 
1.60 
1.80 
2.00 
2.50 
3.00 

4.178 

2.963 

2.400 

2.094 

1.910 

1.791 

1.709 

1.650 

1.560 

1.511 

4.753 

3-517 

2.952 

2.649 

2.468 

2.352 

2.273 

2.217 

2.132 

2.086 

5.406 

4.161 

3.599 

3.301 

3.125 

3.013 

2.937 

2.884 

2.802 

2.759 

6.143 

4.898 

4.343 

4.052 

3.881 

3.773 

3.700 

3.648 

3.569 

3.529 

6.969 

5.732 

5.188 

4.903 

4.737 

4.632 

4.562 

4.512 

4.436 

4.396 

7.817 

6.665 

6.131 

5.853 

5.691 

5.589 

5.518 

5 .'t 71 

5.398 

5.359 


MgSN for cp = 0 

Mo^ 

0 

0.20 

0.40 

0.60 

0.80 

1.00 

0.60 

.80 

1.00 

1.20 

1.40 

1.60 

1.80 

2.00 

2.50 

3.00 

0.5655 

.7274 

.8740 

1.006 

1.123 

1.226 

1.318 

1.397 

1.556 

1.670 

0.5577 

.7070 

.8415 

.9586 

1.060 

l.l 46 

1.220 

1.282 

1.402 

1.483 

0.5432 

.6863 

.8087 

.9120 

.9981 

1.070 

1.129 

1.178 

1.269 

1.328 

0.5319 

.6654 

.7761 

.8665 

.9398 

.9990 

1.047 

1.086 

1.155 

1.199 

0.5203 

.6444 

.7440 

.8229 

.8850 

.9338 

.9727 

1.004 

1.058 

1.091 

0.5095 

.6234 

.7128 

.7815 

.8342 

.8748 

.9064 

.9313 

.9739 

.9996 


(Ms)®=o , 

- — for cp " 0 

\a/ 8 N 

0 

0.20 

0.40 

0.60 

0.80 

1.00 

0.60 

.80 

1.00 

1.20 

i. 4 o 

1.60 

1.80 

2.00 

2.50 

,3.mj 

1.021 

1.021 

1.019 

1.017 

1.016 

1.014 

1.013 

1.012 

1.011 

1.010 

1.024 

1.022 

1.020 

1.018 

1.016 

1.015 

1.014 

1.013 

1.011 

1.010 

1.025 

1.023 

1.020 

1.018 

1.016 

1.015 

1.014 

1.013 

1.012 

1.011 

1.025 

1.023 

1.020 

1.018 

1.016 

1.015 

1.014 

1.013 

1.012 

1.012 

1.026 

1.022 

1.019 

1.017 

1.016 

1.015 

1.014 

1.013 

1.012 

1.012 

1.025 

1.022 

1.019 

1.017 

1.016 

1.015 

1.014 

1.014 

1.013 

1.012 


(Ms)g>=o ^ 

Mn = 30 

\a/ 5 N 

Mq 5 n\ 

0 

0.20 

0.40 

0.60 

0.80 

1.00 

0.60 

.80 

1.00 

1.20 

1.40 

1.60 

1.80 

2.00 

2.50 

3.00 

1.021 

1.021 

1.019 

1.017 

1.016 

1.014 

1.013 

1.012 

1.011 

1.010 

1.021 

1.019 

1.017 

1.014 

1.012 

1.010 

1.008 

1.007 

1.005 

1.004 

1.020 

1.017 

1.013 

1.010 

1.007 

1.005 

1.004 

1.002 

1.000 

.9988 

1.019 

1.014 

1.010 

1.006 

1.004 

1.001 

.9994 

.9981 

.9954 

.9940 

1.017 

1.012 

1.007 

1.003 

.9998 

.9974 

.9955 

.9940 

.9914 

.9897 

1.014 

1.009 

1.004 

.9996 

.9964 

.9939 

.9919 

.9903 

.9875 

.9857 




(Ms), 

Mk 

for cp = 60 

0 


Na/8n 







MoSJK 

0 

0.20 

0.40 

0.60 

0.80 

1.00 

0.60 

1.021 

1.016 

1.009 

1.001 

0.9924 

0.9623 

.80 

1.021 
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Schematic diagram of polar coordinate system for determining supersonic flow 

around o lifting cone. 
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Figure 2- Variation of cone half -angle with free -stream Mach number 
showing range of applicability of cone theory. 
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Figure J— Dimensions of pressure distribution test models showing location of pressure orifices. 
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M a = 0° 




i> = 0° and 180^ 


^ = 45^ and 135° 



^ -- 90° 



— Theory 

o Experiment (measured 
from schHeren 
photographs) 


(b) a =5° 


Figure 4- Shock-wave shapes for an ! 1.42° semivertex angle 
cone at ^ 5.05 and a ^ 0° 5° tO° and 15° 
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Theory 

o Experiment (measured 
from sch/ieren 
photographs) 


(d) a = 15*’ 


Figure 4- Concluded. 
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(a) a -- 




Theory 

o Experiment (measured from 
schHeren photographs) 


(b) a » 5° 


Figure 5- Shock-wave shapes for an 18.92° semivertex angle cone at 
Mo -- 5.05 and a - 0° 5° i0° and !5° 
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=90° 

(c) a = 10 ° 








— Theory 

o Experiment (measured from 
schtieren photographs) 


^ - 90 ° 


(d) a = !5° 


Figure 5 — Concluded. 
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Figure 6- Circumferential variation of pressure coefficient for two test cones at Mo=300 and a = C, and !5? 
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(b) a = 5^ (c) a = io^ (d) a - /5* 

Figure 7.- Circumferential variation of pressure coefficient for two test cones at Mq= 4.25 and a =0^, 5? iOf and 15*: 
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Figure 8- Circumferential variation of pressure coefficient for two test cones at Mo - 5.05 and a = 5f iO't and !5? 
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Figure 9.- Variation of pressure coefficient along ogives at a =0“ and s 3 , 00 , 4.25, and 5.05. 
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Symbol x/i 

o 0 (Cone data) 

° .27 (flagged symbols 

O .54 denote negative 

^ .80 values of 

^ LOO 

Shock- expansion method 

(Stone's second order 
solution employed at 
vertex. (See figure 6)) 
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Meridian angle, </!>, degrees 

(a) Mo = 3.00; K » .60 

Figure iO.- Circumferential variation of pressure coefficient on a fineness ratio 5 ogive at a = 5% IC*, and /5f 
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(b) Mo = 4.25i £5 

Figure !0- Continued, 
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Pressure coefficient, P 




(c) ^ 505i Kr IDI 

Figure !0- Concluded. 
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(a) Mo = 3.00, K=i.00 

Figure it- Circumferential variation of pressure coefficient on a fineness ratio 

3 ogive at a = 5% t01 and /5? 



Pressure coefficient. 
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Symbol x/i 

o O (Cone data) 

□ .22 (flagged symbols 

O .44 denote negative 

^ .66 values of 

^ 1.00 

Shock-expansion method 



Meridian angle, degrees 
(b) Mo = 4.25i K=i.42 


Figure //.- Continued. 
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(c) Mo = 5.05, K = 1.68 


Figure II.- Concluded. 
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Shock- expansion method 

(employing experimentally 
determined cone 
\ conditions) 
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Figure 12.- Circumferential variation of pressure coefficient on a fineness ratio 5 ogive at a = 15°. 
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Figure 13- Circumferential variation of pressure coefficient on a fineness ratio 3 ogive at Of = !5? 
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Figure 14 - Shock-wave shapes for ogives at two values of the similarity parameter, K-, Mq - 5.05, Q = !0^. 
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Figure 15. Pressure distributions predicted by hypersonic slender-body theory for ogives at a =0° 
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Figure 16.— Pressure distributions predicted by hypersonic stender-body theory for ogives at 01=5? 
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(b) K = !.68 


Figure 17- Normal-force coefficients for ogives at Mq = 5.05. 
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